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We study a model in which particles interact with short-ranged attractive and long-ranged re- 
pulsive interactions, in an attempt to model the equilibrium cluster phase recently discovered in 
sterically stabilized colloidal systems in the presence of depletion interactions. At low packing frac- 
tion particles form stable equilibrium clusters which act as building blocks of a cluster fluid. We 
study the possibility that cluster fluids generate a low-density disordered arrested phase, a gel, via 
a glass transition driven by the repulsive interaction. In this model the gel formation is formally 
described with the same physics of the glass formation. 



Gels are disordered arrested systems capable of sus- 
taining weak shear stress. Although commonly used in 
food, cosmetics and detergents, the physical processes 
involved in the transformation of a solution capable of 
flowing into a highly viscous substance, and eventually 
into a gel state, are far from being understood. The gel 
state can be generated via several routes. Interrupted 
phase separation, bond percolation, irreversible diffusion 
limited aggregation, and attractive glass transition have 
been suggested as possible mechanisms for gel formation 
in polymeric and colloidal systems. 

Sterically stabilized colloidal systems have often been 
used as an experimental realization of the hard-sphere 
model, to study both fluid-crystal transition and glass 
formation. The addition of small polymers in solution 
modifies the effective colloid-colloid interaction by deple- 
tion mechanisms, complementing the hard-sphere repul- 
sion with a short-range attraction, on a length scale of 
the order of a few percent of the colloid size. 

At small volume fractions 4> and low temperature T , 
the systems experience dynamic arrest in a transient gel 
structure P, Recent experiments suggest that the 
gel region in the (</> — T) plane is contiguous to a clus- 
ter phase, in which particles H, ^ Hd form freely diffus- 
ing clusters of finite size. Studies of the phase diagram 
of short-ranged (uncharged) attractive colloidal systems 
show that these systems are characterized by fluid and 
crystal phases and by a very flat metastable liquid-gas 
spinodaL which extends to considerably low values of 

0, S Q ■ In simulations of short-ranged attractive po- 
tentials, gel phases have only been observed as a result 
of an arrested phase separation, but with no contiguity 
with an equilibrium cluster phase 

In this Letter we provide evidence, via molecular dy- 
namics (MD) simulations, that when short-range attrac- 
tive interactions are complemented by long-range repul- 
sive interactions, equilibrium cluster phases, as well as 
disordered arrested states contiguous to cluster phases, 
can be generated. The long-range repulsive interactions 
may arise from screened electrostatic interactions or from 



effective interactions accounting for the presence of co- 
solutes in solution piiEli |. 

The possibility of generating a cluster phase, i.e., an 
equilibrium state in which colloidal particles partition 
into stable clusters via the competition between attrac- 
tion (which favors cluster growth at low T) and long- 
range repulsion (which favors low local particle densities 
and, therefore, small aggr egates) has been discussed in 
previous theoretical works ha UJ,hJ. Here we show 
that in the presence of an appropriate balance between 
attraction and repulsion the system forms stable clusters, 
around an optimal cluster size. To quantify this propo- 
sition we calculate the ground state energy of clusters of 
different sizej3 Q for the case of particles interacting 
with a potential V{r) which is the sum of a short-range 
attractive potential Vsr. (modeled for simplicity with the 
generalization to very large a of the Lennard- Jones 2a— a 
potential, as proposed by Vliegenthart et al. j^) 



VsR{r) = 4e 



(1) 



and of a screened electrostatic repulsive interaction Vy 
(modeled by a Yukawa potential) 



Vvir) = A 



(2) 



In the following, the parameters e and a entering in Vsr. 
are chosen as units of energy and distance. The param- 
eter a has been chosen equal to 100, corresponding to 
an attractive range of about three percent of a (see in- 
set in Fig^. For this value of a, the liquid-gas critical 
temperature for Vsr. is found to be 0.235 ± 0.005. Fig. 
shows that the cluster ground state energy, for particles 
interacting with the potential Vsr + Vy for appropri- 
ate values of A and ^ has a minimum at a finite size 
N* , indicating that clusters of size larger than N* are 
energetically disfavored (and, hence, that liquid conden- 
sation is inhibited). Since entropic contributions to the 
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FIG. 1: Ground state energy for connected clusters composed 
of A'^ particles for different choices of A and ^, with a — 100. 
While in the A — case, a monotonous decrease of the energy 
— linear in N~^^^ as expected from the balance between bulk 
and surface effects — is found, for large ^ the lowest energy 
state is reached for a finite A'^* value, which provides an upper 
limit to the size of thermodynamically stable clusters. The 
inset shows Vsr, Vy, and their sum for the case A = 0.05 and 
^ = 2. 



free energy will always favor small clusters, aggregation 
will never proceed beyond N*. The existence of an op- 
timal cluster size (which will of course depend on T and 
(j)) suggests that, at low T, clusters act as building blocks 
of the supra-molecular cluster fluid. 

The cluster-cluster interaction is essentially controlled 
by the Yukawa repulsion, due to the short-range nature of 
the attractive part of the potential. It can be calculated 
analytically, if one approximates the cluster to a sphere of 
radius R composed of a homogeneous density p of points 
interacting via the potential Vy (Eq. Q). The resulting 
cluster-cluster interaction is, again, of the Yukawa type 
with the same ^ but a renormalized amplitude A{R) 



A{R) 



R 
1 



R 



- 1 



(3) 

Therefore, aggregation of particles into clusters of size N 
(and radius R ~ iV^/^) has two effects: i) renormalization 
of the amplitude of the interaction, and ii) decrease of the 
cluster number density n according to 



n(i?) 



1 

N' 



(4) 



Since ^ is not modified, the aggregation process can be 
seen as the generator of a flow, defined via Eqs. (01 and 
Q, in the phase diagram of the pure repulsive Yukawa 
potential. 

Fig. 121 shows an extended phase diagram for the pure 
Yukawa potential. Together with the fluid-crystal coex- 
istence lines, previously calculated by Kremer et al. [Tsj . 
and Meiier and Frenkel [19^ . we show an iso-diffusivity 
curve|23 (i.e., the locus of point with the same diffusion 
constant) and the ideal glass line, evaluated within the 
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FIG. 2: Extended phase diagram for the Yukawa potential, for 
^ = 2, reporting the fiuid-crystal phase boundary lines |l^. 
the iso-diffusivity locus {D = 10~* in unit of ^/A^^Jm) and 
the ideal MCT glass line. The fiow line, calculated according 
to Eqs. Q and describes a system of Yukawa particles 
{A = 0.05, ^ = 2), with packing fraction = Tr/Qa^n = 0.125 
and T — 0.20, assumed to progressively form spherical ho- 
mogeneous clusters of size 1, 2, 4, 8,. . ., 128. When the clus- 
ter size reaches N ~ 47, the melting line is crossed and the 
fluid transforms into a fcc-crystal whose elementary units are 
clusters of that size. If crystallization is avoided, the ideal 
MCT glass line is crossed when the average cluster size reaches 
~ 86. These figures are in very good agreement with the 
simulation data shown in Fig. |3 



Mode Coupling Theory [21| (MCT) using exact numerical 
structure factors as input. The latter curve provides an 
upper limit to the glass transition line. It is important to 
note that the long-range nature of the repulsion enables 
the formation -at low ksT/A values- of solid phases with 
extremely low density (named Wigner crystal or Wigner 
glass HlHil). 

An example of the flow in the Yukawa phase diagram 
induced by clustering is shown in Fig. |21 The starting 
point is assumed to be a system composed of monomers 
[N — 1) interacting via the potential Vsr + Vy- If T is 
sufficiently low, a clustering phenomenon will take place, 
driven by the attractive part of the potential. Such ag- 
gregation process acts as a generator of a flow in the 
phase diagram, moving the system, due to the effective 
increase of the repulsive amplitude (Eq. |2Jl and to the as- 
sociated decrease of the cluster number density (Eq. 0J, 
closer and closer to the melting and glass lines. If the 
optimal size is reached before crossing the melting line, 
the equilibrium state for the system will correspond to a 
fluid phase of clusters. In the other case, an arrest phe- 
nomenon will take place, driven by the crystallization or 
-more likely because of the inherent self-generated cluster 
polydispersity- by a glass transition. With this mecha- 
nism arrested disordered states with small can be gen- 



3 




FIG. 3: Cluster size distribution, on a log-log scale, for (j> — 
0.125. For clarity, curves for different T have been shifted 
along the vertical axis progressively by one decade. The inset 
shows the T-dependence of the average cluster size. It also 
shows, for each T, the size at which the crystallization and 
the ideal glass transition lines are crossed, according to the 
flow described in Fig. |21 



erated. This new route to the gel state demonstrates, in 
the case of weakly charged short-ranged attractive col- 
loids, the identity of the gel and the glass states of mat- 
ter. 

We next provide numerical evidence that the above 
proposed mechanism does lead to disordered arrested 
states contiguous to a cluster phase. We perform MD 
simulations of particles interacting with the potential 
VsR + Vy 01 , both along a constant T and along a con- 
stant (j) path. Here, we focus on the specific choice of 
A = 0.05 and ^ = 2, referring to future publications for 
a more extensive analysis. This choice of parameters is 
characteristic of colloidal particles in weakly polar sol- 
vents Hillll. 

The cluster size distribution, shown in Fig. |31 pro- 
vides a characterization of the clustering process. While 
at high T the cluster size is monotonically decreasing, 
around ksT/e = 0.2 the distributions develop a peak at 
a finite value, signaling the presence of stable equilib- 
rium aggregates composed of 20 — 30 particles. We note 
that such resulting cluster phase can be considered as 
an echo of the interrupted spinodal decomposition, since 
clustering takes places just below the liquid-gas critical 
temperature of the bare Vsr potential. 

For each of the studied T, the unclustered system 
{N = 1) can be represented as a point in Fig. |21 A 
flow curve, starting from each of these points, can be 
calculated according to Eqs.^land^ In this way, a the- 
oretical prediction for the average cluster size at which 
crystallization and (MCT) vitrification are expected can 
be derived. Such curves are shown in the inset together 
with the T-dependence of the average cluster size calcu- 
lated from the MD data. For the present choice of A, ^ 
and (f), structural arrest (either in the form of crystalliza- 
tion or in the form of vitrification) is predicted when the 
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FIG. 4: Evolution of static and dynamic quantities on cooling 
at (/) = 0.125. (Top) S'(g). Note the development of a peak 
at small q, reflecting the formation of clusters in the system. 
(Bottom) F{q,t) at qa = 2.7. At T = 0.1 the autocorrelation 
function does not decay to zero any longer, highlighting the 
formation of a structurally arrested phase. 



average cluster size is > 35 particles {ksT/e < 0.15). 

Before turning to the evaluation of dynamical prop- 
erties, we discuss the evolution of the structure with T, 
analyzing the static structure factor S{q). As shown in 
Figure EJa, S{q) progressively develops a peak at wave- 
vectors q associated with distances in real space of the 
order of several particle diameters, reflecting the near- 
est neighbor cluster-cluster distance. The amplitude of 
this novel peak changes significantly in a small tempera- 
ture range around the T-region where the fiow diagram 
predicts the crossing of the crystallization line. There- 
fore, the system changes from a solution of interacting 
particles toward a system of larger and larger interacting 
clusters. In the cluster phase, the location of the peak 
changes very weakly both with T and with 0. It is inter- 
esting to note that both the location and the amplitude 
of the cluster-cluster peak are in agreement with the ex- 
perimental data for colloidal gels reported in Ref. . To 
confirm that the aggregation process indeed generates a 
disordered arrested state, we show in Fig. 01 (Bottom) the 
coherent scattering function F{q,t). It displays a (one- 
step) decay, with a characteristic time which increases on 
cooling. At the lowest studied T, the correlation func- 
tion does not decay to zero, signahng that the system is 
trapped in a non-ergodic state. The transition from fluid 
to glass is very sharp, again in full agreement with results 
for colloidal gels Q . 

In summary, we have proposed and studied a model 
which shows equilibrium cluster phases and arrested 
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states at low densities, opening the way for accurate 
comparisons between experiments and theoretical pre- 
dictions. Contrary to what has been proposed before 
for attractive colloids, we suggest that, in these systems, 
the gel phase is not stabilized by the short-range attrac- 
tion (and hence its formation is not connected to the 
attractive glass transition observed at high packing frac- 
tions ^J). The simulation data (Figs.|3|and21) show that 
the competition between isotropic short-range attraction 
and long-range repulsion — of the ty pe f ound in colloidal 
systems and in protein solutions 29', 3(| — is capable of 
generating a very rich range of phenomena |3l| . Within 
the present model, the formation of low density arrested 



states can be modeled as a glass transition, where clusters 
(as opposed to particles) are trapped in cages generated 
by the long-range repulsions. It will be interesting to 
study if the arrested disordered state is metastable with 
respect to a a fee or bcc crystal of "spherical" clusters (a 
Wigner crystal) or with respect to lamellar phases, which 
has been observed in related models [13, [ij, . 
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tractive colloidal systems. Authors acknowledge support 
from MIUR Colin 2002 and Firb. 
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